Structure of the human KMN complex and implications for regulation of its assembly.

Biorientation of chromosomes during cell division is necessary for precise dispatching of a mother cell's chromosomes into its two daughters. Kinetochores, large layered structures built on specialized chromosome loci named centromeres, promote biorientation by binding and sensing spindle microtubules. One of the outer layer main components is a ten-subunit assembly comprising Knl1C, Mis12C and Ndc80C (KMN) subcomplexes. The KMN is highly elongated and docks on kinetochores and microtubules through interfaces at its opposite extremes. Here, we combine cryogenic electron microscopy reconstructions and AlphaFold2 predictions to generate a model of the human KMN that reveals all intra-KMN interfaces. We identify and functionally validate two interaction interfaces that link Mis12C to Ndc80C and Knl1C. Through targeted interference experiments, we demonstrate that this mutual organization strongly stabilizes the KMN assembly. Our work thus reports a comprehensive structural and functional analysis of this part of the kinetochore microtubule-binding machinery and elucidates the path of connections from the chromatin-bound components to the force-generating components.

Substituting density functional theory in reaction barrier calculations for hydrogen atom transfer in proteins.

Hydrogen atom transfer (HAT) reactions are important in many biological systems. As these reactions are hard to observe experimentally, it is of high interest to shed light on them using simulations. Here, we present a machine learning model based on graph neural networks for the prediction of energy barriers of HAT reactions in proteins. As input, the model uses exclusively non-optimized structures as obtained from classical simulations. It was trained on more than 17 000 energy barriers calculated using hybrid density functional theory. We built and evaluated the model in the context of HAT in collagen, but we show that the same workflow can easily be applied to HAT reactions in other biological or synthetic polymers. We obtain for relevant reactions (small reaction distances) a model with good predictive power (R2 ∼ 0.9 and mean absolute error of <3 kcal mol-1). As the inference speed is high, this model enables evaluations of dozens of chemical situations within seconds. When combined with molecular dynamics in a kinetic Monte-Carlo scheme, the model paves the way toward reactive simulations.

Mechanical forces control the valency of the malaria adhesin VAR2CSA by exposing cryptic glycan binding sites.

Plasmodium falciparum (Pf) is responsible for the most lethal form of malaria. VAR2CSA is an adhesin protein expressed by this parasite at the membrane of infected erythrocytes for attachment to the placenta, leading to pregnancy-associated malaria. VAR2CSA is a large 355 kDa multidomain protein composed of nine extracellular domains, a transmembrane helix, and an intracellular domain. VAR2CSA binds to Chondroitin Sulphate A (CSA) of the proteoglycan matrix of the placenta. Shear flow, as the one occurring in blood, has been shown to enhance the (VAR2CSA-mediated) adhesion of Pf-infected erythrocytes on the CSA-matrix. However, the underlying molecular mechanism governing this enhancement has remained elusive. Here, we address this question by using equilibrium, force-probe, and docking-based molecular dynamics simulations. We subjected the VAR2CSA protein-CSA sugar complex to a force mimicking the tensile force exerted on this system due to the shear of the flowing blood. We show that upon this force exertion, VAR2CSA undergoes a large opening conformational transition before the CSA sugar chain dissociates from its main binding site. This preferential order of events is caused by the orientation of the molecule during elongation, as well as the strong electrostatic attraction of the sugar to the main protein binding site. Upon opening, two additional cryptic CSA binding sites get exposed and a functional dodecameric CSA molecule can be stably accommodated at these force-exposed positions. Thus, our results suggest that mechanical forces increase the avidity of VAR2CSA by turning it from a monovalent to a multivalent state. We propose this to be the molecular cause of the observed shear-enhanced adherence. Mechanical control of the valency of VAR2CSA is an intriguing hypothesis that can be tested experimentally and which is of relevance for the understanding of the malaria infection and for the development of anti placental-malaria vaccines targeting VAR2CSA.

The role of flow in the self-assembly of dragline spider silk proteins.

Hydrodynamic flow in the spider duct induces conformational changes in dragline spider silk proteins (spidroins) and drives their assembly, but the underlying physical mechanisms are still elusive. Here we address this challenging multiscale problem with a complementary strategy of atomistic and coarse-grained molecular dynamics simulations with uniform flow. The conformational changes at the molecular level were analyzed for single-tethered spider silk peptides. Uniform flow leads to coiled-to-stretch transitions and pushes alanine residues into ß sheet and poly-proline II conformations. Coarse-grained simulations of the assembly process of multiple semi-flexible block copolymers using multi-particle collision dynamics reveal that the spidroins aggregate faster but into low-order assemblies when they are less extended. At medium-to-large peptide extensions (50%-80%), assembly slows down and becomes reversible with frequent association and dissociation events, whereas spidroin alignment increases and alanine repeats form ordered regions. Our work highlights the role of flow in guiding silk self-assembly into tough fibers by enhancing alignment and kinetic reversibility, a mechanism likely relevant also for other proteins whose function depends on hydrodynamic flow.

Myristoyl's dual role in allosterically regulating and localizing Abl kinase.

c-Abl kinase, a key signaling hub in many biological processes ranging from cell development to proliferation, is tightly regulated by two inhibitory Src homology domains. An N-terminal myristoyl modification can bind to a hydrophobic pocket in the kinase C-lobe, which stabilizes the autoinhibitory assembly. Activation is triggered by myristoyl release. We used molecular dynamics simulations to show how both myristoyl and the Src homology domains are required to impose the full inhibitory effect on the kinase domain and reveal the allosteric transmission pathway at residue-level resolution. Importantly, we find myristoyl insertion into a membrane to thermodynamically compete with binding to c-Abl. Myristoyl thus not only localizes the protein to the cellular membrane, but membrane attachment at the same time enhances activation of c-Abl by stabilizing its preactivated state. Our data put forward a model in which lipidation tightly couples kinase localization and regulation, a scheme that currently appears to be unique for this non-receptor tyrosine kinase.

Allosteric activation of vinculin by talin.

The talin-vinculin axis is a key mechanosensing component of cellular focal adhesions. How talin and vinculin respond to forces and regulate one another remains unclear. By combining single-molecule magnetic tweezers experiments, Molecular Dynamics simulations, actin-bundling assays, and adhesion assembly experiments in live cells, we here describe a two-ways allosteric network within vinculin as a regulator of the talin-vinculin interaction. We directly observe a maturation process of vinculin upon talin binding, which reinforces the binding to talin at a rate of 0.03 s-1. This allosteric transition can compete with force-induced dissociation of vinculin from talin only at forces up to 10 pN. Mimicking the allosteric activation by mutation yields a vinculin molecule that bundles actin and localizes to focal adhesions in a force-independent manner. Hence, the allosteric switch confines talin-vinculin interactions and focal adhesion build-up to intermediate force levels. The 'allosteric vinculin mutant' is a valuable molecular tool to further dissect the mechanical and biochemical signalling circuits at focal adhesions and elsewhere.

DOPA Residues Endow Collagen with Radical Scavenging Capacity.

Here we uncover collagen, the main structural protein of all connective tissues, as a redox-active material. We identify dihydroxyphenylalanine (DOPA) residues, post-translational oxidation products of tyrosine residues, to be common in collagen derived from different connective tissues. We observe that these DOPA residues endow collagen with substantial radical scavenging capacity. When reducing radicals, DOPA residues work as redox relay: they convert to the quinone and generate hydrogen peroxide. In this dual function, DOPA outcompetes its amino acid precursors and ascorbic acid. Our results establish DOPA residues as redox-active side chains of collagens, probably protecting connective tissues against radicals formed under mechanical stress and/or inflammation.

Collagen breaks at weak sacrificial bonds taming its mechanoradicals.

Collagen is a force-bearing, hierarchical structural protein important to all connective tissue. In tendon collagen, high load even below macroscopic failure level creates mechanoradicals by homolytic bond scission, similar to polymers. The location and type of initial rupture sites critically decide on both the mechanical and chemical impact of these micro-ruptures on the tissue, but are yet to be explored. We here use scale-bridging simulations supported by gel electrophoresis and mass spectrometry to determine breakage points in collagen. We find collagen crosslinks, as opposed to the backbone, to harbor the weakest bonds, with one particular bond in trivalent crosslinks as the most dominant rupture site. We identify this bond as sacrificial, rupturing prior to other bonds while maintaining the material's integrity. Also, collagen's weak bonds funnel ruptures such that the potentially harmful mechanoradicals are readily stabilized. Our results suggest this unique failure mode of collagen to be tailored towards combatting an early onset of macroscopic failure and material ageing.

Disulfide bond reduction and exchange in C4 domain of von Willebrand factor undermines platelet binding.

BACKGROUND: The von Willebrand factor (VWF) is a key player in regulating hemostasis through adhesion of platelets to sites of vascular injury. It is a large, multi-domain, mechano-sensitive protein that is stabilized by a net of disulfide bridges. Binding to platelet integrin is achieved by the VWF-C4 domain, which exhibits a fixed fold, even under conditions of severe mechanical stress, but only if critical internal disulfide bonds are closed. OBJECTIVE: To determine the oxidation state of disulfide bridges in the C4 domain of VWF and implications for VWF's platelet binding function. METHODS: We combined classical molecular dynamics and quantum mechanical simulations, mass spectrometry, site-directed mutagenesis, and platelet binding assays. RESULTS: We show that 2 disulfide bonds in the VWF-C4 domain, namely the 2 major force-bearing ones, are partially reduced in human blood. Reduction leads to pronounced conformational changes within C4 that considerably affect the accessibility of the integrin-binding motif, and thereby impair integrin-mediated platelet binding. We also reveal that reduced species in the C4 domain undergo specific thiol/disulfide exchanges with the remaining disulfide bridges, in a process in which mechanical force may increase the proximity of specific reactant cysteines, further trapping C4 in a state of low integrin-binding propensity. We identify a multitude of redox states in all 6 VWF-C domains, suggesting disulfide bond reduction and swapping to be a general theme. CONCLUSIONS: Our data suggests a mechanism in which disulfide bonds dynamically swap cysteine partners and control the interaction of VWF with integrin and potentially other partners, thereby critically influencing its hemostatic function.

Mechano-Redox Control of Mac-1 De-Adhesion by PDI Promotes Directional Movement Under Flow.

BACKGROUND: Neutrophil migration is critical to the initiation and resolution of inflammation. Macrophage-1 antigen (Mac-1; CD11b/CD18, αMß2) is a leukocyte integrin essential for firm adhesion to endothelial ICAM-1 (intercellular adhesion molecule 1) and migration of neutrophils in the shear forces of the circulation. PDI (protein disulfide isomerase) has been reported to influence neutrophil adhesion and migration. We aimed to elucidate the molecular mechanism of PDI control of Mac-1 affinity for ICAM-1 during neutrophil migration under fluid shear. METHODS: Neutrophils isolated from whole blood were perfused over microfluidic chips coated with ICAM-1. Colocalization of Mac-1 and PDI on neutrophils was visualized by fluorescently labeled antibodies and confocal microscopy. The redox state of Mac-1 disulfide bonds was mapped by differential cysteine alkylation and mass spectrometry. Wild-type or disulfide mutant Mac-1 was expressed recombinantly in Baby Hamster Kidney cells to measure ligand affinity. Mac-1 conformations were measured by conformation-specific antibodies and molecular dynamics simulations. Neutrophils crawling on immobilized ICAM-1 were measured in presence of oxidized or reduced PDI, and the effect of PDI inhibition using isoquercetin on neutrophil crawling on inflamed endothelial cells was examined. Migration indices in the X- and Y-direction were determined and the crawling speed was calculated. RESULTS: PDI colocalized with high-affinity Mac-1 at the trailing edge of stimulated neutrophils when crawling on ICAM-1 under fluid shear. PDI cleaved 2 allosteric disulfide bonds, C169-C176 and C224-C264, in the ßI domain of the ß2 subunit, and cleavage of the C224-C264 disulfide bond selectively controls Mac-1 disengagement from ICAM-1 under fluid shear. Molecular dynamics simulations and conformation-specific antibodies reveal that cleavage of the C224-C264 bond induces conformational change and mechanical stress in the ßI domain. This allosterically alters the exposure of an αI domain epitope associated with a shift of Mac-1 to a lower-affinity state. These molecular events promote neutrophil motility in the direction of flow at high shear stress. Inhibition of PDI by isoquercetin reduces neutrophil migration in the direction of flow on endothelial cells during inflammation. CONCLUSIONS: Shear-dependent PDI cleavage of the neutrophil Mac-1 C224-C264 disulfide bond triggers Mac-1 de-adherence from ICAM-1 at the trailing edge of the cell and enables directional movement of neutrophils during inflammation.

FERM domains recruit ample PI(4,5)P2s to form extensive protein-membrane attachments.

The four-point-one ezrin-radixin-moesin homology (FERM) protein domain is a multifunctional protein-lipid binding site, constituting an integral part of numerous membrane-associated proteins. Its interaction with the lipid phosphatidylinositol-4,5-bisphosphate (PIP2), located at the inner leaflet of eukaryotic plasma membranes, is important for localization, anchorage, and activation of FERM-containing proteins. FERM-PIP2 complexes structurally determined so far exclusively feature a 1:1 binding stoichiometry of protein and lipid, with a few basic FERM residues neutralizing the -4 charge of the bound PIP2. Whether this picture from static crystal structures also applies to the dynamic interaction of FERM domains on PIP2 membranes is unknown. We here quantified the stoichiometry of FERM-PIP2 binding in a lipid bilayer using atomistic molecular dynamics simulations and experiments on solid supported membranes for the FERM domains of focal adhesion kinase and ezrin. In contrast to the structural data, we find much higher average stoichiometries of FERM-PIP2 binding, amounting to 1:3 or 1:4 ratios, respectively. In simulations, the full set of basic residues at the membrane interface, 7 and 15 residues for focal adhesion kinase and ezrin, respectively, engages in PIP2 interactions. In addition, Na ions enter the FERM-membrane binding interface, compensating negative PIP2 charges in case of high charge surpluses from bound PIP2. We propose the multivalent binding of FERM domains to PIP2 in lipid bilayers to significantly enhance the stability of FERM-membrane binding and to render the FERM-membrane linkage highly adjustable.

Intrinsically disordered region of talin's FERM domain functions as an initial PIP2 recognition site.

Focal adhesions (FAs) mediate the interaction of the cytoskeleton with the extracellular matrix in a highly dynamic fashion. Talin is a central regulator, adaptor protein, and mechano-sensor of FA complexes. For recruitment and firm attachment at FAs, talin's N-terminal FERM domain binds to phosphatidylinositol 4,5-bisphosphate (PIP2)-enriched membranes. A newly published autoinhibitory structure of talin-1, where the known PIP2 interaction sites are covered up, lead us to hypothesize that a hitherto less examined loop insertion of the FERM domain acts as an additional and initial site of contact. We evaluated direct interactions of talin-1 with a PIP2 membrane by means of atomistic molecular dynamics simulations. We show that this unstructured, 33-residue-long loop strongly interacts with PIP2 and can facilitate further membrane contacts, including the canonical PIP2 interactions, by serving as a flexible membrane anchor. Under force as present at FAs, the extensible FERM loop ensures talin maintains membrane contacts when pulled away from the membrane by up to 7 nm. We identify key basic residues of the anchor mediating the highly dynamic talin-membrane interaction. Our results put forward an intrinsically disordered loop as a key and highly adaptable PIP2 recognition site of talin and potentially other PIP2-binding mechano-proteins.

Hydropersulfides inhibit lipid peroxidation and ferroptosis by scavenging radicals.

Ferroptosis is a type of cell death caused by radical-driven lipid peroxidation, leading to membrane damage and rupture. Here we show that enzymatically produced sulfane sulfur (S0) species, specifically hydropersulfides, scavenge endogenously generated free radicals and, thereby, suppress lipid peroxidation and ferroptosis. By providing sulfur for S0 biosynthesis, cysteine can support ferroptosis resistance independently of the canonical GPX4 pathway. Our results further suggest that hydropersulfides terminate radical chain reactions through the formation and self-recombination of perthiyl radicals. The autocatalytic regeneration of hydropersulfides may explain why low micromolar concentrations of persulfides suffice to produce potent cytoprotective effects on a background of millimolar concentrations of glutathione. We propose that increased S0 biosynthesis is an adaptive cellular response to radical-driven lipid peroxidation, potentially representing a primordial radical protection system.

Bond dissociation energies of X-H bonds in proteins.

Knowledge of reliable X-H bond dissociation energies (X = C, N, O, S) for amino acids in proteins is key for studying the radical chemistry of proteins. X-H bond dissociation energies of model dipeptides were computed using the isodesmic reaction method at the BMK/6-31+G(2df,p) and G4(MP2)-6X levels of theory. The density functional theory values agree well with the composite-level calculations. By this high level of theory, combined with a careful choice of reference compounds and peptide model systems, our work provides a highly valuable data set of bond dissociation energies with unprecedented accuracy and comprehensiveness. It will likely prove useful to predict protein biochemistry involving radicals, e.g., by machine learning.

Structure and dynamics of the von Willebrand Factor C6 domain.

Von Willebrand disease (VWD) is a bleeding disorder with different levels of severity. VWD-associated mutations are located in the von Willebrand factor (VWF) gene, coding for the large multidomain plasma protein VWF with essential roles in hemostasis and thrombosis. On the one hand, a variety of mutations in the C-domains of VWF are associated with increased bleeding upon vascular injury. On the other hand, VWF gain-of-function (GOF) mutations in the C4 domain have recently been identified, which induce an increased risk of myocardial infarction. Mechanistic insights into how these mutations affect the molecular behavior of VWF are scarce and holistic approaches are challenging due to the multidomain and multimeric character of this large protein. Here, we determine the structure and dynamics of the C6 domain and the single nucleotide polymorphism (SNP) variant G2705R in C6 by combining nuclear magnetic resonance spectroscopy, molecular dynamics simulations and aggregometry. Our findings indicate that this mutation mostly destabilizes VWF by leading to a more pronounced hinging between both subdomains of C6. Hemostatic parameters of variant G2705R are close to normal under static conditions, but the missense mutation results in a gain-of-function under flow conditions, due to decreased VWF stem stability. Together with the fact that two C4 variants also exhibit GOF characteristics, our data underline the importance of the VWF stem region in VWF's hemostatic activity and the risk of mutation-associated prothrombotic properties in VWF C-domain variants due to altered stem dynamics.

A conformational transition of the D'D3 domain primes von Willebrand factor for multimerization.

Von Willebrand factor (VWF) is a multimeric plasma glycoprotein that is critically involved in hemostasis. Biosynthesis of long VWF concatemers in the endoplasmic reticulum and the trans-Golgi is still not fully understood. We use the single-molecule force spectroscopy technique magnetic tweezers to analyze a previously hypothesized conformational change in the D'D3 domain crucial for VWF multimerization. We find that the interface formed by submodules C8-3, TIL3, and E3 wrapping around VWD3 can open and expose 2 buried cysteines, Cys1099 and Cys1142, that are vital for multimerization. By characterizing the conformational change at varying levels of force, we can quantify the kinetics of the transition and stability of the interface. We find a pronounced destabilization of the interface on lowering the pH from 7.4 to 6.2 and 5.5. This is consistent with initiation of the conformational change that enables VWF multimerization at the D'D3 domain by a decrease in pH in the trans-Golgi network and Weibel-Palade bodies. Furthermore, we find a stabilization of the interface in the presence of coagulation factor VIII, providing evidence for a previously hypothesized binding site in submodule C8-3. Our findings highlight the critical role of the D'D3 domain in VWF biosynthesis and function, and we anticipate our methodology to be applicable to study other, similar conformational changes in VWF and beyond.

Martini 3 coarse-grained force field for poly(para-phenylene ethynylene)s.

Poly(para-phenylene ethynylene)s, or short PPEs, are a class of conjugated and semi-flexible polymers with a strongly delocalized π electron system and increased chain stiffness. Due to this, PPEs have a wide range of technological applications. Although the material properties of single-chains or mixtures of few PPE chains have been studied in detail, the properties of large assemblies remain to be fully explored. Here, we developed a coarse-grained model for PPEs with the Martini 3 force field to enable computational studies of PPEs in large-scale assembly. We used an optimization geometrical approach to take the shape of the π conjugated backbone into account and also applied an additional angular potential to tune the mechanical bending stiffness of the polymer. Our Martini 3 model reproduces key structural and thermodynamic observables of single PPE chains and mixtures, such as persistence length, density, packing and stacking. We show that chain entanglement increases with the expense of nematic ordering with growing PPE chain length. With the Martini 3 PPE model at hand, we are now able to cover large spatio-temporal scales and thereby to uncover key aspects for the structural organization of PPE bulk systems. The model is also predicted to be of high applicability to investigate out-of-equilibrium behavior of PPEs under mechanical force.

ATP allosterically stabilizes integrin-linked kinase for efficient force generation.

SignificanceThe pseudokinase integrin-linked kinase (ILK) is a central component of focal adhesions, cytoplasmic multiprotein complexes that integrate and transduce biochemical and mechanical signals from the extracellular environment into the cell and vice versa. However, the precise molecular functions, particularly the mechanosensory properties of ILK and the significance of retained adenosine triphosphate (ATP) binding, are still unclear. Combining molecular-dynamics simulations with cell biology, we establish a role for ATP binding to pseudokinases. We find that ATP promotes the structural stability of ILK, allosterically influences the interaction between ILK and its binding partner parvin at adhesions, and enhances the mechanoresistance of this complex. On the cellular level, ATP binding facilitates efficient traction force buildup, focal adhesion stabilization, and efficient cell migration.

Mechanical force can enhance c-Src kinase activity by impairing autoinhibition.

Cellular mechanosensing is pivotal for virtually all biological processes, and many molecular mechano-sensors and their way of function are being uncovered. In this work, we suggest that c-Src kinase acts as a direct mechano-sensor. c-Src is responsible for, among others, cell proliferation, and shows increased activity in stretched cells. In its native state, c-Src has little basal activity, because its kinase domain binds to an SH2 and SH3 domain. However, it is known that c-Src can bind to p130Cas, through which force can be transmitted to the membrane. Using molecular dynamics simulations, we show that force acting between the membrane-bound N-terminus of the SH3 domain and p130Cas induces partial SH3 unfolding, thereby impeding rebinding of the kinase domain onto SH2/SH3 and effectively enhancing kinase activity. Forces involved in this process are slightly lower or similar to the forces required to pull out c-Src from the membrane through the myristoyl linker, and key interactions involved in this anchoring are salt bridges between negative lipids and nearby basic residues in c-Src. Thus, c-Src appears to be a candidate for an intriguing mechanosensing mechanism of impaired kinase inhibition, which can be potentially tuned by membrane composition and other environmental factors.

Accurate Free Energies for Complex Condensed-Phase Reactions Using an Artificial Neural Network Corrected DFTB/MM Methodology.

Semiempirical methods like density functional tight-binding (DFTB) allow extensive phase space sampling, making it possible to generate free energy surfaces of complex reactions in condensed-phase environments. Such a high efficiency often comes at the cost of reduced accuracy, which may be improved by developing a specific reaction parametrization (SRP) for the particular molecular system. Thiol-disulfide exchange is a nucleophilic substitution reaction that occurs in a large class of proteins. Its proper description requires a high-level ab initio method, while DFT-GAA and hybrid functionals were shown to be inadequate, and so is DFTB due to its DFT-GGA descent. We develop an SRP for thiol-disulfide exchange based on an artificial neural network (ANN) implementation in the DFTB+ software and compare its performance to that of a standard SRP approach applied to DFTB. As an application, we use both new DFTB-SRP as components of a QM/MM scheme to investigate thiol-disulfide exchange in two molecular complexes: a solvated model system and a blood protein. Demonstrating the strengths of the methodology, highly accurate free energy surfaces are generated at a low cost, as the augmentation of DFTB with an ANN only adds a small computational overhead.